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ABSTRACT 

Aims. Substellar objects have extremely long life-spans. The cosmological consequence for older objects are low abundances of heavy 
elements, which results in a wide distribution of objects over metallicity, hence over age. Within their cool atmosphere, dust clouds 
become a dominant feature, affecting the opacity and the remaining gas phase abundance of heavy elements. We investigate the influ- 
ence of the stellar metallicity on the dust formation in substellar atmospheres and on the dust cloud structure and its feedback on the 
atmosphere. This work has implications for the general question of star formation and of dust formation in the early universe. 
Methods. We utilize numerical simulations in which we solve a set of moment equations in order to determine the quasi-static dust 
cloud structure (Drift). These equations model the nucleation, the kinetic growth of composite particles, their evaporation and the 
gravitational settling as a stationary dust formation process. Element conservation equations augment this system of equations includ- 
ing the element replenishment by convective overshooting. The integration with an atmosphere code (Phoenix) allows to determine 
a consistent {T, p, v^onv)- structure {T - local temperature, p - local pressure, v^onv - convective velocity), and, hence, also to calculate 
synthetic spectra. 

Results. A grid of Drift-Phoenix model atmospheres was calculated for a wide range of metallicity, [M/H] 6 [-1-0.5,-0.0,-0.5,.. .,-6.0], 
to allow for a systematic study of atmospheric cloud structures throughout the evolution of the universe. We find dust clouds in even 
the most metal-poor ([M/H] =-6.0) atmosphere of brown dwarfs. Only the most massive among the youngest brown dwarfs and giant 
gas planets can resist dust formation. For very low heavy element abundances, a temperature inversion develops which has a drastic 
impact on the dust cloud structure. 

Conclusions. The combination of metal depletion by dust formation and the uncertainty of interior element abundances makes the 
modeling of substellar atmospheres an intricate problem in particular for old substellar objects. We further show that the dust-to-gas 
ratio does not scale linearly with the object's [M/H] for a given effective temperature. The mean grain sizes and the composition of 
the grains change depending on [M/H] which influences the dust opacity that determines radiative heating and cooling as well as the 
spectral appearance. 

Key words, astrochemistry - Methods: numerical - Stars: atmospheres - Stars: low-mass, brown dwarfs 



!>■ 1. Introduction limit for stellar formati on is stil l not set tled. For solar abun- 

dances, simulations by IWhitwo rth & St amatellosI (i2006) pre- 
>j Condensation of supersaturated gas species becomes a major is- diet a limit for the mass of coll apsing gas c lou ds as l ow as 
^ sue in ultracool atmospheres, i.e. that of brown dwarfs and giant 0.001 to 0.004Mo. Results by e.g. lBossI dioOlh a ndlBatd (l2005l) 



planets. The formation of solid dust particles or liquid droplets fall in a sir nilar mas s regim e. Observations by iGreaves et all 

depletes the gas phase of elements and, hence, strongly affects (|2003h and iGreavesI (l2005h present evidence for the forma- 

the chemical composition of the remaining gas phase. In addi- tion of objects with masses below the deuterium burning 

tion, the forming dust particles represent an important opacity Umit via a gravitational collapsing molecular clou d fragment. 

source, which ultimately results in a significant backwarming The p r esence of discs arou n d S-lOMjup-objects dNatta etal] 

and an inc reased gas temper ature compared to a cloud-free at- ( |2002l) . iLuhman et all (|2005|) . IScholz & Javawardhanal (l2008l) ) 

mosphere jTsuii et alJll996al) . is also evidence for the formation of such low-mass objects 

Besides the dust formation induced element deficiency, the via gravitational clou d collaps. According to Omukai (200j) 

substellar object can already have formed from a metal-poor gas, and IClark et alj (l2008l) subsolar mass objects may already have 

hence, being intrinsically metal-deficient. The first generation of formed at [M/H]>-6.0, though heating effects may have ob- 

stars was comprised of extremely massive objects. These stars structed the collapse of low mass gas cl ouds and depopulated 

were responsible for a strong metal enrichment of the interstel- the range between [M/H]=-4.0 and -5.0. iJappsen et all (l2009h 

lar medium, from which new stars formed. The coiTesponding note that the transition of the early universe IMF to the present 

enrichment of the ISM with more and more heavy elements had day IMF at around [M/H] = 10 - 10"^ seem to be determined 

an important impact on the star formation, as it brought about the by dust cooling which in turn is determined by the grain size 

collap se of much smaller cloud fragments (Smith & Sigurdsson distribution and composition of the dust forming at these low 



collapsi 
|2007[) . 



Hence, succeeding star generations contained 'growing metalicities. Those extreme metallicities can be attributed to a 

numbers of smaller objects. The debate about the lower mass much earlier state of the universe, when the continuously run- 
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ning heavy element enrichment of the interstellar medium by 
nucleosynthesis inside massive stars and supernovae had not 
progressed as far as today. The implication is a wide distribu- 
tion of objects over the metallicity range. However, all but the 
lowest mass objects of high age, i.e. of extremely low metal- 
licities, have akeady exceeded their life-span as main-sequence 
stars. Therefore, only ultracool stars and substellar objects re- 
main from those earlier stages of the universe. Hence, substellar 
spectra need to be characterised not only by effective temper- 
ature, Teff, and surface gravity, log(g), but also by the interior, 
well-mixed element abundances as relic of the conditions at their 
formation. 

The presence of dust in ultracool dwarf and substella r at- 
mospheres has been suggested first by iLunine et al.l (Il986l) and 
was noticed by iJones & Tsuiil (Il997h . First observational evi- 
dence for c onvec tive a ctivity in substellar objects was found by 
iNoU et all (1 19971) and Oppenhei mer et al.l (119981) . By then, the 
first b rown dwarfs had jus t been discovered by iRebolo et al.l 
( Il995h and iNakaiima etaTI (Il995h . The obs ervation of meta l- 
poor ultracool stars in the galactic halo fe.g.. lChristliebl (l2006l) ) 
led to t he classifica tion into dwarfs, subdwarfs and extreme sub- 
dwarfs (|G izis"1997). A later revision of this system was done by 
iLepine et al., (.2007,) . introducing the even more metal-poor ultra- 
su bdwarfs class. The fir st subste llar subdwarfs were discovered 
bv lBurgasser et al.1 (l2003i) and Burgasseij (|2004|) . 

As the importance of dust clouds in late-M-, L- and 
T-type dwarf atmospheres had been recognised, models 
were developed to help to understand these objects. Initially, 
the clouds were parameterised by a constant grain size 
dTsuji et al.l ll996b: Burrows et al. 2006), later time-scale argu- 
ments were applied (Allardetal. 2001), or grai n size distri- 
butions guided by Earth-observa tions were used (ICooper et al.l 
120031: 1 Ackerman & Marlevll200ll) . All these models apply com- 
plete (S = 1) or almost (S = 1.01) phase-equilibrium. A de- 
tailed compa rison study of current dust models was made by 
iHelUng et al.l (2Q08a). 

We consider a stationary dust formation process, where seeds 
form high up in the atmosphere from a highly supersaturated 
gas, grow to macroscopic particles of around a jjm in size, gravi- 
tationally settle into deeper layers and eventually evaporate as 
the local temperature be comes too high for thermal stability 
dWoitke & Hellindl2003l) . The basic idea is that a "dirty" soUd 
mantle will grow on top of the seed particles. This dirty mantle 
is assumed to be composed of numerous small islands of dif- 
ferent pure condensates. The formatio n of islands is sup ported 
by experiments in solid state physics dLedieu et al. |2005|) and 
by observations of coated t errestrial dust particles dLevin et al.l 
[19961 iKorhonen et al.ll2003h . 



2. Method 

We make use o f the Drift-Ph oenix model atmosph e re code , 
introduced by iDehnI (|2007|) and iHelUng et al.] d2008bl) . 

The general-purpos e mo d el atmosph ere code Phoenix 

dHauschildt & BaronI 119991: iBaron et al.1 l2003h solves the 
gas-phase equation of state, provides the atmosphere structure 
to Drift, determines the gas o pacities and solves the r adiative 
ti-ansfer. The Drift code by iHelling. Woitke. & Thil (l2008cl) 
is included as a module in order to calculate the dust clouds, 
which feed back on both the thermodynamic al structures and 
the radiation field. An iteration of this method allows the deter- 
mination of the atmosphere and dust cloud properties and yields 
the corresponding synthetic spectra. For a given atmosphere 
structure (T, p, v'conv), provided by Phoenix, the phase non- 



equilibrium dust formation with subsequent precipitation and 
element replenishment by convective overshooting is calculated 
by Drift. The opacities of the composite dust particles, which 
are required for the radiati ve transfer, are det ermined by using 
eff'ective medium theorv (|Bruggemanl 1935b and Mie theory 
dMielll908t IWolf & Voshchinnikovll2004 . All dust species ai-e 
included as opacity source since our dust grains are made of a 
mixtures of these compounts. 



2. 1 . Gas-phase chemistry 

The chemistry in Phoenix considers the 40 most important el- 
ements in several ionization states and 47 molecule species. 
Though larger chemi cal systems are possible, the results dif- 
fer only marg inally dDehnll2007l) . For the solar element abun- 
dances, defining [M/H 1-0.0 in our models, we utilize data from 
iGrevesse et alj dl992h . We are aware that the elemental com- 
position of brown dwarfs and planets is dependent on the ini- 
tial composition of the gas cloud they formed from, which it- 
self depends critically on the internal physics of its local pre- 
ceding star generations, like, e.g., the efficiency of the various 
burning cycles and the resulting amount of elements. Especially 
for extremely low metallicities, the element abundance pattern is 
completely different to the one we know from the solar enviro n- 
ment (e.g., iFrebel et alj d2005h . ICampbell & Lattanziol d2008')). 
However, we neglect variations of the individual element abun- 
dances relative to each other and their impact on metallicity 
model sequences for the sake of simplicity and determine the 
varied metallicity values by a single scaling-factor. These values 
are varied according to the metallicity under consideration such 
that [M/H]=-h0.5 ...-6.0. 

The Phoenix chemistry utilizes the altitude-dependent gas 
phase abundances ei of the elements Mg, Si, Ti, O, Fe, and 
Al provided by Drift, i.e., it takes into account those element 
abundances which have been altered throughout the atmosphere 
by depletion (nucleation, dust growth), enrichment (evapora- 
tion) and redistribution (precipitation, convective overshooting). 
In the DRiFT-module, we calculate the number densities of all 
gaseous species, including the number density of the key reac- 
tant n^*^^ as described in IWoitke & HellingI d2004l) according to 
pressure, temperature and the calculated depth-dependent ele- 
ment abundances e, in chemical equilibrium. For the well-mixed, 
deep element abundances e!' we use the scaled solar abundances 
of the model input. For those elements that are not included in 
the Drift chemistry calculations, we set e, = e'^. 



3. Growth reactions 

3.1. Dust seed formation 

The nucleation ra te Jj, is calculated for (Ti0 2)Af-clusters accord- 
ing to Eq. (34) in ' Helling & Woi tke (2006"), applying the modi- 
fied classical nucleation theory of Gail et al. ( 1984). We use the 
value of th e surface tension cr fitted to small cluster da ta by 
iJeong et al.l dIOOQ) as outlined in lWoitke & Helling! (|2004|) . 

3.2. Dust growth/evaporation & refractive indices 

Th e dust growth of dirty p articles is modelle d according 
to 'Woitke & Helling diool . iHelUng & Woife e ('200 ^ and 
Helling et al. (2008c). However, compared to Hellin get al.l 
d2008cl) . we consider only the 7 most important solids (Ti02[s], 
Al203[s], Fe[s], Si02[s], MgO[s], MgSi03[s], Mg2Si04[s]) 
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Solid s 


Surface reaction r 


Key species 


TiO,[s] 


TiO, — I Ti02[s] 


Ti02 


rulile 


Ti + 2 H,0 — . Ti02[s] + 2 H, 


Ti 


(1) 


TiO + H,0^ Ti0,[s] + H, 


TiO 


(A) 


TiS + 2 H,0 Ti6,[s] + HjS + H, 


TiS 


SiOiM 


SiO, SiO,[s] 


Si02 


silica 


SiO + H2O — > Si02[s] + H, 


SiO 


(3)(B) 


SiS + 2 H2O —> Si6,[s] + H2S + H2 


SiS 


Fe[s] 


Fc Fe[s] 


Fe 


solid iron 


FcO + H, — > Fe[s] + H2O 


FeO 


(1) 


FcS + H, ^ Fe[s] + HjS 


FeS 


(B) 


Fc(0H)2 + H2 ^ Fc[s] + 2 H2O 


Fe(OH)2 


MgO[s] 
periclase 

(3) 

(B) 


MgO MgO[s] 
Mg + H,0 — 1 MgO[s] + H, 
2 MgOH — > 2 MgO[s] + H, 
Mg(0H)2 — > Mgb[s] + H2O 


MgO 
Mg 

i;,MgOH 
Mg(OH)2 


MgSi03[s] 


Mg + SiO + 2 H,0 MgSiO,[s] + H, 


minlMg, SiO| 


enstatite 


Mg + SiS + 3 H,0 MgSiO,[s] + H,S + 2 H, 


minlMg, SiS| 


(3) 
(C) 


2 MgOH + 2 SiO + 2 H,0 2 MgSi6,[s] + 3 H, 

2 MgOH + 2 SiS + 2 H2O — > 2 MgSiO,[s] + 2 H,S + 2 H, 

Mg(OH), + SiO ^ 2 MgSiO,[s] + H, 

Mg(0H)2 + SiS + H2O MgSi03[s]"+ H2S+ H2 


minlV-MgOH, i/,SiO| 
minlV-MgOH, i/,SiS| 
min|Mg(0H)2, SiO| 
min(Mg(0H)2, SiS| 


Mg,Si04[s] 


2 Mg + SiO + 3 H2O — > Mg2Si04[s] + 3 H2 


min( 'iMg, SiO| 


forsterite 
(3) 
(C) 


2 MgOH + SiO + H2O Mg2Si04[s] + 2 H, 
2 Mg(0H)2 + SiO Mg2Si04[s] + H2O + H2 
2 Mg + SiS + 4 H2O — > Mg2Si04[s] + H2S + 3 H2 


min|>/,MgOH, SiO| 
min|"tMg(OH)2, SiO| 
miniv-Mg, SiS] 




2 MgOH + SiS + 2 H2O — rMg2Si04[s] + H2S + 2 H2 
2 Mg(0H)2 + SiS — > Mg2Si04[s] + H2 + H2S 


min|'/,MgOH. SiS| 
min('tMg(0H)2, SiS| 


AhOjM 


2 AlOH + H2O Al20,[s] + 2 H2 


|/,A10H 


aluminia 


2 AIH + 3 H2O Al,0,[s] + 4 H, 


'/,A1H 


(3) 


A1,0 + 2 H2O — > Al,b,[s] + 2 H2 


A1,0 


(D) 


2 AIS + 3 H,0 — 1 AUO,[s] + 2 H,S + H, 


|/,A1S 




2 AIO,H — > Al,0,[s] + H,0 


|/,A10,H 



Table 1. Chemical surface reactions r assumed to form the 
solid materials s. The efficiency of the reaction is limited by the 
collision rate of the key species. The notation 'k in the r.h.s. col- 
umn means that only every second collision (and sticking) event 
initiates one reaction. Data sources fo r the supersaturation ratio s 
(an d saturation vapou r pressures): (l)'Hellin2 & Woitke' ('20061); 
(2) lNuth & FergusonI (12006); (3) Sharp & Huebner ( 1990). Data 
sou rces for the pure s olid refractive indices: (A) Posch (2008); 
(B')|Pang(fT98l|l99l: (O .Jageret aL(.2003.) : (D) .Palika99il) . 
iBegemann et all(ll997l) . 



made of 6 different elements for which we solve the result- 
ing, stiff sys t em of 15 (3x Eq.(8), 7x Eq.(9), 5xEq.(10) in 
iHelUng et al.l (l2008d) ^ dust moment and element conservation 
equations. 

The 7 considered solids are formed by 32 chemical surface 
reaction (see Table [T]i in phase-non-equilibrium to calculate the 
formation and composition of the dirty grains. Our selection is 
guided by the most stable condensates which yet have simple 
stoichiometric ratios, ensuring that these solids can be easily 
built up from the gas phase (f or more details see Sect. 3.4. in 
iHell ing. Woit ke. & Thil(l2008cl) ). Note that we consider dust for- 
mation by gas-solid reactions only and omit solid-solid reactions 
and lattice rearrangements i nside the grains in our model (see 
iHelUng & RietmeiieJ (120091) regarding lattice rearrangements in 
clouds in substellar atmospheres). Table[T]also contains the solid 
refractive index data used in the Mie calculations. 



3.3. Element replenishment 

A truly static atmos phere would not contain any dust 
dWoitke & Helling||2004i) . Therefore, we include mixing by con- 
vection and overshooting by assuming an exponential decrease 
of th e mass exchange frequency in the radiative zone (Eq. 9 in 
dmitke & Helling 2004) with yS = 2.2 and = 2/(//pVconv)), 
which serves to replenish the upper atmospheres and keeps the 
cycle of dust formation running. 
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Fig. 1. (r, /7)-structures for three model sequences of vary- 
ing metalHcity [M/H]: Top: Teff=2000K, log(g)=5.0; Center: 
Teft=2500K, log(g)=5.0; Buttom: Teif=2000K, log(g)=3.0 



3.4. Grain size distribution function 

The dust opacity calculations with effective medium theory 
and Mie theory require the grain size distribution function 
/(fl) [cm"'*] at every depth in the atmosphere, where a [cm] is 
the grain radius. This function is not a direct result of the dust 
moment method. Only the total dust particle number density 
rid - pLo and the mean particle size {a) = y/T/iATrj Li / Lq are 
direct results that have been used in Helling & Woitke (200^). 
In this paper, we reconstruct f(a) from the calculated dust 
moments Lj (j = 1 ...4) in an approxima t e way as described 
in Appendix A in iHelling. Woitke. & Thil (l2008cl) . In this pa- 
per, we apply the double delta-peaked size distribution function 
/(fl) = A^i 6{a -ai)+N2 6(a - 02) where 6 is the Dirac-function, 
Ni,N2 [cm"-'] are two dust particle densities and 01,02 [cm] are 
the two corresponding particle radii. 
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4. Results: The model grid 

In this section we will describe the influence of metallicity on 
the model atmospheres and the dust cloud. This will be done 
for three showcase model sequences to exemplarily cover the 
parameter space for effective temperature and surface gravity 
which covers brown dwarfs and giant gas planets. We think that 
such a detailed demonstration of our dust model results is nec- 
essary as dust cloud model for brown dwarfs differ considerably 
with respect to dust properties (Helling et al..,2008ai) . 
The sequences are calculated for: 

Teff = 2000/: and log(g)=5.0, 
Teff = 2500/: and log(g)=5.0 
Teff = 2000/: and log(g)=3.0, 

covering the metalUcity values [M/H] e [-1-0.5, -0.0,-0.5,...,-6.0]. 
In the following, we will concentrate on only [M/H]e[-i-0. 5,-0.0,- 
1.5,-4.0,-6.0]. 

4.1. Metallicity-dependent temperature-pressure structures 

At first, we will take a look how the temperature-pressure struc- 
tures are affected by the metallicity (Fig.[T]i. As one might ex- 
pect, the local temperature decreases with [M/H] for a fixed 
pressure, due to the smaller gas opacities. However, the outer 
atmospheres of the low [M/H] models feature a strong tempera- 
ture inversion. It is caused by a tilt in the chemical equilibrium, 
favouring the formation of methane instead of carbon monoxide, 
which frees up oxygen for the formation of water and metal ox- 
ides (e.g., Allard & Hauschildt (1995)). The inversion becomes 
stronger for higher effective temperatures and higher surface 
gravities for a given [M/H] -value. 

Figure[T] demonstrates that the local pressure increases for a 
given local temperature with decreasing metallicity but with in- 
creasing gravity. This relative pressure increase strengthens the 
contribution of strongly pressure-dependent gas-phase opacity 
sources like the collision induced absorption but also the con- 
tribution of the pressure-broadened alkali lines. iFortnev et al.l 
(1200 8) present similar findings for cooler substellar objects with 
Teff ^ 1400K. Hence, the effect of increasing pressure with de- 
creasing metal content holds across the whole substellar range of 
Teff. Therefore, clouds in giant planets but also in low-metallicity 
substellar objects form at even lower gas temperatures at a com- 
parable pressure (see Sect. 14. 3t . 

The dust clouds represent a strong opacity source. Therefore, 
the backwarming caused by the dust particles strongly increases 
the local gas temperature at and below the cloud location. Hence, 
dense cloud layers coincide with a strong inward temperature 
gradient. Especially the high [M/H] models feature two bends 
around 1500K and 2000K. The drop of the local tempera- 
ture gradient is caused by the starting evaporation of the solid 
species. The lower temperature bend is attributed to the evap- 
oration of the silicate species (MgO[s], Si02[s], MgSi03[s], 
Mg2Si04[s]), while the second corresponds to evaporation of the 
thermally more stable species (Fe[s], Ti02[s], Al203[s]). 

In the Teff=2000K model for log(g)=5.0 and [M/H]=-4.0, 
there is a bump in the (T,p)-structure around 6-10^ dyncm"^. 
It akeady appears for [M/H] =-2.0 and remains until [M/H]=- 
5.0. The origin is the superposition of two effects, namely the 
decreasing deep interior abundances and the decreasing deple- 
tion of the gas phase with the metallicity. This causes a locally 
higher absolute Ti abundance in the gas phase compared to at- 
mospheres of even higher metallicities. The opacity by the re- 



sulting increased amount of TiO is the reason for this feature in 
the (T,p)-structure. 

4.2. Metallicity influence on the (vconv,/')-sfrucfure 

Solar-metallicity models ([M/H]=0.0) suggest that the dust 
cloud location and the convectively unstable zone are well sep- 
arated. Even so, convection provides a mechanism for element 
replenishment of the upper atmosphere layer, which assures sta- 
tionarity of the cloud formation process. The so-called convec- 
tive o vershooting is paramete rised in Drift by a mixing-time 
scale (IWoitke & Helling|l2004 following the idea of momentum 
conservation of the convective mass elements (see lLudwig et al.l 
(l2002i |2006|) ) whose motion does not instantly end where the 
Schwarzschild criterion is no longer met, but decays slowly into 
the higher layers ( Woitke & Helli ng 2004). A higher maximum 
convection velocity consequently yields a more efficient mixing 
into the upper layers. 

Due to the vanishing local opacities with decreasing [M/H], 
the radiative flux is less obstructed and, hence, the convective en- 
ergy transport becomes less important. For this reason, the con- 
vection zone is shifted inwards to higher pressures and densities 
for decreasing metallicities. According to mixing length theory, 
the resulting mean convection velocity Vconv decreases, which 
results in a less efficient convective overshooting. Because of a 
weaker radiation field for decreasing effective temperatures and 
higher gas densities for increasing surface gravity, one finds a 
similar effect on the convective velocity and overshooting for 
these parameters. 

Hence, the convective overshooting is more efficient in hot, 
low-gravity objects. An improved mixing consequently enables 
the dust cloud to persist into higher atmosphere layers. Hence, 
hot-Jupi ters should have hig h-altitu de cloud s which supports 
what Richardson et al.l (l2007i) and Pont et al.l (12008 ) suggested 
for their observed giant gas planets. 

4.3. Metallicity-dependent cloud structure 

We study the influence of the stellar metallicity on the dust cloud 
structure to provide insight into possible diversities of cloud 
properties for substellar objects originating from different evolu- 
tionary states of the universe, i.e. for objects which have formed 
from interstellar gas clouds of different metal-enrichment. 

Figure|2] contains the results including the (T, /?)-structures 
for better comparison. Depicted are the number of dust particles 
forming per volume element and time, i.e. the nucleation rate 
7, [cm^^s"'] (2"'' panel), the net gr owth ve l ocit y = YjsX^ 
[cms-'] (3"* pan el; Eq. (7) in HeUin g et all (|20b8c)), the mean 
grain size {a} = -^fTfiAn) L\ I L() [cm] (4* panel), and the equilib- 
rium drift velocity (vdi) = ■\/7rgpd{a} I i^pcr) [cms"'] (5* panel), 
which is determined for gravitationally settling particles of mean 
size, retarded by friction for the case of free molecular flows 
(Woitke & Helling.2003.) . 

4.3.1. The general, high-[M/H] cloud structure 

As is shown in Fig.|2] all high metallicity models exhibit a typ- 
ical cloud structure which was already shown in our previous 
publications. At the highest altitudes w ith its dust-poor depleted 
gas (compare 'Woi tke & HelUn j (|2004|) ). dust seeds form which 
are unable to grow perceptibly due to the extremely low densi- 
ties at these altitudes. The mean grain size, {a}, remains of the 
order of the seed particle size. Due to the low density, the num- 
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10° 10^ lO* 10° 10° 10° 10' 10' 10^ 10* 10' 10° 10' 10"' 10"' 10° 10' 10' 10' 10* 10' 

P,os [dyn cm-'] 

[M/H] = + 0.5 [M/H] = -0.0 [M/H] = -1.5 [M/H] = -4.0 [M/H] = -6.0 

Fig. 2. Cloud structures for three model sequences of varying metallicity [M/H]. Left: Teff=2000K, log(g)=5.0; Center: 
Teff=2500K, log(g)=5.0; Right: Teff=2000K, log(g)=3.0. 1"' row: local temperature T [K], l""^ row: nucleation rate 7, [cm^^s^'], 
'y^ row: net growth velocity net [cms"'], 4* row: mean grain size {a) [cm], 5* row: mean drift velocity (vdr) [cms"']. 



ber of the dust particles remains very low (not shown), but the 
particles can settle very fast into deeper layer for the same rea- 
son. The net grain growth velocity, Xnei, increases proportionally 
to the gas density, while the drift velocity, (vdr), decreases. From 
a certain depth downward, the mean particle size increases from 
the initial size of the seeds. Together with the increasing nucle- 
ation rate, the dust growth starts to cause a steep drop of the local 
gas-phase element abundances with depth (compare Fig.|6]). The 
consequence is that the growth velocity eventually decreases to 
a degree, that the particles settle too fast to sustain a percep- 
tible growth of the mean grain radius at this low-pressure re- 
gion. Gravitational settling does not influence the process of seed 
formation since nucleation proceeds much faster. However, the 
strong gas-phase depletion will impede the formation of more 
seeds as less material is available. As the precipitating parti- 
cles drop into atmosphere layers which are closer to the edge 
of the convection zone, the convective upmixing becomes more 
efficient. The more efficient mixing results in a continuously in- 
creasing nucleation rate with depth as more material is available 
to form seed from. The gas-phase depletion by the exponential 
increase of the number of new seed particles compensates for 
the on-going growth of the grains. Hence, the mean grain radius 
remains constant, even though there is growth, which is strong 
enough to cause a exponentially decreasing element abundances, 
shown in Fig.|6l In this region, the dust cloud becomes optically 
thick enough to cause a perceptible backwarming of up to several 
hundred degrees (Fig.[TJ. As the local temperatures are increas- 
ing and the abundances are decreasing with depth, the nucleation 
rate reaches its apex and drops to zero very fast over the next few 
atmospheric layers. This also means, the particle production has 
reached its maximum. Hence, dust existing in deeper layers must 
have rained in from above. Therefore, changes in the number 
density of the dust particles in these deeper atmospheric layers 



are caused by accumulation, drift and evaporation. Interestingly, 
the dust number density remains of the same order compared to 
the site of the nucleation maximum, but the grain size distribu- 
tion changes dramatically. 

The end of nucleation sparks a strong increase of the mean 
grain size. The frictional force caused by the high densities at 
these layers has strongly decreased the drift velocity. Thus, it 
takes the dust particles much longer to cross a given atmospheric 
layer. This allows the particles to capture more and more reac- 
tants from the gas phase, resulting in an even stronger growth. 
In turn, the mean grain radius causes an increase of the drift 
velocity. This is balanced by the increasing gas density, which 
is why we observe a plateau-like feature in the drift velocity 
(Fig. 121 lowest panel). The numerous yum-sized particles below 
the nucleation maximum cause a strong backwarming. 

As soon as the evaporation of a certain species exceeds its 
growth, it is vanishing from the solid phase. This is accompanied 
by a decreasing and possibly even negative growth velocity. That 
is the reason for the observed kinks (Fig.|2]l in the mean grain 
radius curves and all consequential kinks in other quantities like, 
e.g., the effective temperature, as mention in Sec. l4.1l As soon as 
the last solid species starts to evaporate, the dust cloud typically 
vanishes very fast with depth. 

4.3.2. The influence of metallicity on tine cloud structure 

The convective mixing becomes less efficient with decreasing 
[M/H] (Sec. |4.2| l. Thus, the dust clouds are located at higher 
gas pressures and densities for lower [M/H]. However, this is 
also caused by an additional factor, which is a temperature in- 
version in the outer atmosphere (Sec. 14. 11 Fig. [1]), developing 
for metallicities below [M/H]=-3.0. Within this inversion zone, 
the gas phase is no longer supersaturated for the lowest consid- 



6 



S. Witte et al.: Dust in brown dwarfs and extra-solar planets 



ered metallicities. Thus, no dust can form there and the typical 
structure of undisturbed dust clouds is lost. In these model at- 
mospheres, the dust cloud is bottled up in a temperature "valley" 
between the inversion zone and the hotter inner atmosphere. This 
valley becomes more destinctive for increasing surface gravity, 
i.e., it will vanish for planet-like objects but persist within the 
coolest amongst the brown dwarfs (Fig.[Tli. 

Depending on metallicity, the atmospheres can have one 
or two zones of efficient nucleation, i.e. two nucleation rate 
maxima 7™"". High-metallicity models show only the outer 
nucleation zone, as rep orted by Woitke & Helling ( 2003') and 
iHelhng. Woitke. & Thil (l2008cl) . Very metal-deficient models 
([M/H]<-4.0) develop a very shallow outer nucleation zone and 
an additional inner zone of much more efficient seed forma- 
tion. The outer nucleation zone is already indicative for a layer 
of haze. Nucleation takes place throughout the dust cloud for 
low metallicities in contrast to the higher metallicity models, 
where nucleation ceases in the middle of the cloud. Generally, 
the maximum nucleation rate is decreasing with metallicity 
down to [M/H]=-3.0, even though the maximum is shifted in- 
wards to higher densities. The sole exception to this trend is 
the [M/H] =-1-0.5 atmosphere model in the hottest sequence pre- 
sented here (Teff=2500K), which counterintuitively shows a 
smaller than the lower-metallicity models of that sequence. 
Since the gas temperature is akeady high in this model, not much 
backwarming is needed to move the model out of the tempera- 
ture window of efficient nucleation (compare and 2"'* panel 
in Fig|2]i. 

For lower metallicities ([M/H]< -3.0), the trend of decreas- 
ing nucleation maximum with [M/H] is reversed due to the ap- 
pearance of the second, inner 7™^". The reason for this is the 
strongly decreasing local temperature above the convection zone 
with decreasing metallicity, which allows a strong supersatura- 
tion (compare Sect. |4.4t . The result is a growing second max- 
imum of the nucleation rate with decreasing metallicity at the 
cloud base. 

Our results show that dust may persist down to extremly low 
values of metallicity. A combination of low metallicity, high ef- 
fective temperature and low surface gravity is required to prevent 
the formation of dust clouds. Therefore, only the log(g)=3.0 se- 
quence, representing the planetary atmospheres in our sample, 
features no dust at all for [M/H]<-4.0. 

The growth velocity of the mean particle typically has two 
local maxima. The outer maximum of the growth velocity de- 
creases and is shifted to higher pressures for decreasing metal- 
licities. It vanishes as soon as the temperature inversion sets in 
for [M/H]< -4.0. The deeper located second maximum of the 
growth velocity is also shifted inwards. However, it is neither 
affected by the temperature inversion, nor does it decrease con- 
tinuously with the metallicity, because the inward shift brings it 
much closer to the convection zone and therefore a more efficient 
element replenishment. Around [M/H] =-4.0 the cloud indeed 
dips into the convection zone for our young giant planet-like 
models (Teft=2000K, log(g)=3.0), which sparks an additional 
growth by one order of magnitude in those models. 

Due to the proximity to the convection zone and the result- 
ing high growth velocity, the maximum mean grain radius in- 
creases slightly with decreasing metallicity down to [M/H]=-3.5. 
Only for lower metallicities, the decreasing abundances are fi- 
nally able to cause a decreasing maximum mean grain size. 

The drift velocity decreases much stronger for lower effec- 
tive temperatures and to a lesser extend for higher surface grav- 
ities. As a lower drift velocity results in a stronger accumula- 



tion of dust particles, the number density in the lower cloud is 
strongly increasing with decreasing effective temperature. 

4.4. The phase-non-equilibrium of the clouds 

The thermodynamic conditions in an undepleted substellar at- 
mosphere guarantee a strong supersaturation of a wide range of 
solid compounds. An initial high supersaturation does not au- 
tomatically yield a phase equilibrium between the depleted at- 
mosphere and the dust cloud as such jHelling. Woitke. & TTiil 
l2008cl) . The phase equilibrium would require a sufficiently high 
number density of reactants as well as high collision and reaction 
rates. However, this is usually not the case for most of the solids 
of interest. Therefore, all three quantities need to be considered 
in order to avoid the overrated dust growth and gas phase deple- 
tion of the phase eq uilibrium a pproach. According to ou r model 
(see IWoitke& Hel ling (20031): iHelling & Woitker(l2006h ). most 
of the atmosphere is in phase non-equilibrium, a chemical equiv- 
alent to the NLTE state of radiation fields. 

Solid grains or liquid droples are composed of so-called 
monomer units. Experiments show that most of the complex 
monomers like Mg2Si0 4 do not exist in the gas phase (e.g. 
iRietmeiier et alj (Il999b ). The formation of such solid units 
would require multi body reactions which makes the probabil- 
ity for their formation in the gas neglegibly small. On the other 
hand, if the reactants are physically absorbed by a preexisting 
surface, their limited phase space permits a much more effi- 
cient growth reaction. The reactant of the least physically ab- 
sorbed amount is called the key species. Its absorption rate de- 
termines the rate of the considered chemical surface reaction 
and, he nce, also the su persaturation ratio for each surface re- 
action jHelling & Woitke 2006). In order to discuss the satu- 
ration of the atmosphere, we apply the effective supersatura- 
tion ratio fo r each considered solid compound as defined in 
Helling, Woit ke. & Thil (l2008cl) . 

The effective supersaturation ratios for each of the consid- 
ered solid species for the three model sequences are shown in 
Fig. [3] For the higher metallicities, the outer atmosphere has 
extremely high supersaturation ratios. The ratio is increasing 
with the pressure, because the concentration of reactants in- 
creases with the pressure. Further inwards, the supersaturation 
ratio reaches a maximum where the dust growth starts to cause 
a strong depletion of the gas phase. Below that point, the su- 
persaturation decreases proportionally to the dwindling element 
abundances. At a c ertain layer, each of the ratios r each es unity. 
In agr eement with iHelhng. Woitke. & Thil (l2008cl) and iHeUingl 
(12009') we find that all solids which contain high-abundant el- 
ements (Mg, Si, O, Fe) achieve phase-equilibrium in the pres- 
sure interval where the cloud has it's maximum dust-to-gas ratio, 
Pd/Pgas (compare Figs. |3] and |5]l, before the supersaturation ratios 
drop further below unity. Solids which contain rare elements (Ti, 
Al) do not feature a phase-equilibrium region. 

For decreasing abundances, a higher gas density is required 
to reach the same amount of supersaturation. However, the su- 
persaturation becomes weaker for lower metallicities. This, su- 
perimposed with the temperature inversion zone of the low 
metallicity models, results in a fast drop below unity of all su- 
persaturation ratios, which means no dust is able to form there. 
Fortunate for the dust cloud in the low metallicity models is the 
decreasing gas temperature between the inversion zone and the 
convection zone. Even though the element abundances are de- 
creasing with the metallicity, the low temperatures bring about 
a larger concentration of reactants. Therefore, the supersatura- 
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tion ratio in the lower atmosphere is increasing with decreasing 
metallicity. 

In the [M/H]=-6.0 model of the log(g)=3.0 sequence, 
Ti02[s] is never sufficiently supersaturated for nucleation to take 
place. Therefore, this model is dust free. 



4.5. Metallicity-dependence of the cloud composition 

Pure silicate dust grains are mostly transparent. Ti02[s] and 
Al203[s] also do not represent strong opacities. However, even 
a small impu rity can turn the transparency into opaqueness 
(IWoitkell2006 h. This means, there is a huge difference between 
a mixture of pure grains and mixture of "dirty" grains. Hence, 
together with the grain size distribution, the composition of the 
mean particle is of special importance for the radiative transfer 
and, thus, for the spectrum. 

Drift was designed to treat the formation of dirty grains. 
Therefore, it considers grains consisting of islands of mixed 
species. The composition changes with atmospheric depth and 
the stellar parameters. The results for the three examples se- 
quences are shown in Fig.|4] As Ti02[s] is the solid species 
which forms the dust seeds, the particles in the outer cloud are 
solely made of it. Ti is rare compared to the other five considered 
elements. Thus, as soon as growth becomes an issue, the volume 
fraction of Ti02[s] decreases fast. At the higher cloud layers, 
where the temperatures are lower, the silicates with traces of iron 
dominate the dust grains. The temperatures increase with depth, 
eventually leading to the evaporation of the silicates, while the 
growth of iron becomes much more efficient, resulting in very 
opaque particles. Slightly deeper, Al203[s] growth dominates 
the volume of the average grains. 

We find that the inward shift of the dust cloud yields a 
stronger maximum of the iron volume fraction. The opposite 
happens with the Al203[s] maximum, which becomes weaker 



The chemical equilibrium for the lower metallicities is 
changed, which leads to a significant preference of MgO[s] over 
the other siHcates for [M/H]<-4.0 (Fig.gl left column). This ef- 
fect is most pronounced in objects with strong dust formation 
and the accompanying strong gas phase depletion. For higher 
effective temperatures (Fig.|4] middle column), it can not be ob- 
served in our models due to the much weaker depletion, which 
yields a delay of the chemical equilibrium shift to even lower 
metallicities. In the lower temperature, low gravity sequence 
(Fig. El right column) the effect is simply not observable because 
no dust is present in the respective low metallicity models for 
the given effective temperature of 2000K. However, the changed 
chemical equilibrium can be observed for these low gravities as 
soon as the temperatures permit the formation of dust. Hence, 
the change from silicate to MgO[s] dominated dust grains for ex- 
tremely low metallicities is observable for cooler objects. Higher 
surface gravities emphasize the extend of the shift. 

The appearing temperature inversion zone in the low metal- 
licity models for high gravities causes the disappearance of the 
wide silicate dominated cloud region, leaving an almost pure 
Fe[s] and Al203[s] cloud. 



4.6. The atmospheric dust content 

The dust to gas ratio of the mass densities, pd/Pgas, is a good 
measure for the amount of dust within the atmosphere (Fig.lSj. 
A lower metallicity typically yields a lower pd/Pgas- The excep- 
tion to this is again the Tefr=2500K [M/H] =-1-0.5 model, because 
of the high temperatures (Sec. |4.3t . One might assume that the 
dust content would scale linearly with the element abundances. 
However, it is not that simple, because the observed inward shift 
of the dust clouds with decreasing metallicity yields a higher 
concentration of condensable elements. This density effect par- 
tially counteracts the decreasing abundances due to decreasing 
metallicities. It becomes more pronounced for cooler, low grav- 
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ity objects. For example, the maximum dust content decreases 
by merely 2.5 orders of magnitude between [M/H] =+0.5 and 
[M/H]=-4.0 in the Teff =2000K log(g)=5.0 sequence. Only when 
the cloud is not able to sink any deeper this trend ends and the 
dust content is finally approaching a linear scaling with the abun- 
dances. 

4. 7. Remaining gas pliase abundances 

Dust formation reduces the local element abundances in the gas 
phase, until nucleation and growth become so weak that the de- 
pletion is fully balanced by the replenishment. At its minimum, 
the gas phase abundances e, of the more rare elements can reach 
values below 10 of their deep interior abundance e, q. Titanium 
is the least abundant of the six presently considered condens- 
ing elements. Especially in the outer atmosphere, which is only 
weakly affected by convective overshooting but its conditions 
allow seed formation, titanium is extremely depleted (Fig.|6]l. 

As the efficiency of dust growth increases inwards, the de- 
pletion gets stronger. Eventually, this ends as soon as the in- 
wards increasing element consumption by the dust formation 
stops to outbalance the increase of the gas density. This is due to 
the inwards decreasing supersaturation of the gas and the con- 
sequently less efficient growth. For supersaturation ratios below 
unity, the respective solids evaporate, returning their constituents 
to the gas phase. This enrichment can reach values of twice the 
deep interior model abundances. 

The depletion follows the shifts of the dust cloud with the 
stellar parameters. There are two opposing influences on the dust 
growth, which cause changes in the depletion. A deeper located 
cloud features a more efficient growth because of the higher den- 
sity. However, the inwards shift caused by a decreasing metallic- 
ity is accompanied by a decreasing growth efficiency. Hence, 
there is a metallicity for which the latter effect starts to can- 
cel the former, resulting in a turn from increasing to decreasing 
maximum value of the relative depletion with decreasing metal- 



licity. For Teff=2000K and log(g) =5.0 we find this maximum 
depletion around [M/H] =-1.0. For higher effective temperature 
the turning point metallicity value is increasing, while it is de- 
creasing for lower surface gravity. 

We conclude that the strong depletion, observed in our mod- 
els, can have a tremendous influence on absorption line profiles. 
For wavelengths at which the dust cloud is mostly transparent, 
one is able to see the deep and undepleted atmosphere layers. 
Due to their high density these layers will dominate the respec- 
tive spectral lines. Because these layers are not depleted by dust 
formation, the full element abundances will be seen in the spec- 
tral lines. However, for the wavelengths at which the dust cloud 
is opaque, the observed abundances will not resemble the well- 
mixed interior abundances of the object. Depending on the wave- 
length, the observed abundances can be several orders of magni- 
tude lower than the actual ones. Therefore, it is extremely com- 
plex to derive individual stellar abundances from spectral lines 
in dust-forming environments in general. 

4.8. Dust content and grain size above t=1 

In the previous sections, the spatial distribution of dust within 
substellar atmospheres was analysed. It remains to be verified 
how much dust is observationally accessible, e.g., by polarime- 
try. In Fig.|7]we show the local gas pressure at the maximum of 
the dust-to-gas ratio P(max{pd/Pg)), the corresponding optical 
depth T and mean grain size {a} for all three model sequences 
(compare to Fig.|2]and|5]l. 

Figure|2] shows optical depth values for one single wave- 
length sample (1.2yum; ti 2;ini) and for the Planck-, Rosseland- 
, mean intensity- and surface-normal flux-weighted opacities 
{tb, Tioss, tj and rp, respectively). At 1.2/im the atmospheres 
are highly transparent. In contrast, the weighted optical depth 
means, in particular tb which in case of the considered cool 
atmospheres emphasizes on the near-infrared where the dust is 
most opaque, are typically larger by one to two orders of magni- 



S. Witte et al.: Dust in brown dwarfs and extra-solar planets 



9 




-2 
-4 


-2 
-4 



: -2 
-4 


-2 
-4 


-2 

-4 
0.05 
0.00 

-0.05 
-0.10 




T,||=2000K, log(g)=5.0 





T„|=2500K, log(g)=5.0 



10' 



10' 1 0-' 1 0" 

[M/H] = + 0.5 



10= 



10' 10' 

[M/H] = 



. Fe 






-0.0 



10' 10' 1 o' 
P505 [dyn cm-'] 

[M/H] = -1.5 



10= 




-2 
-4 


-2 
-4 


-2 
-4 


-2 
-4 


-2 

-4 
0.05 
0.00 

-0.05 
-0.10 



= 2000K, log(g) = 5.0 





Fe 


J,, / ; J - 




T[ 






SI ^ 


I •'- ' 








Al 






Mg ^ 


/ / * * 

•si;, / » 




- X. 


; ; ■ 









10"' 10"' 10° 10' 10' 10' 10" 10= 

[M/H] = -4.0 [M/H] = -6.0 



Fig. 6. Element abundances, £[, for three model sequences of varying metallicity [M/H]. Left: Teff=2000K, log(g)=5.0; Center: 
Teff=2500K, log(g)=5.0; Right: Teff=2000K, log(g)=3.0 



tude. The Rosseland mean is the lowest of the considered mean 
optical depths, as it puts emphasis on the wavelengths of lowest 
opacity. 

The dust-to-gas ratio maximum, max{pd/Pg), is shifted in- 
wards for decreasing metallicities down to [M/H]=-3.0. . .-4.0, 
with the exact value depending on the stellar parameters. 
For lower values, the pd/Pg -maximum rises in altitude, again 
(Fig-Ell- Therefore, all curves in Fig.|2] feature a sharp turn in 
local pressure around these metallicities. This clearly demon- 
states the non-linear influence of the metallicity on the dust con- 
tent in substellar atmospheres. Between about [M/H]=-0.5 and 
+ 1.0 the mean grain size and the optical depths remain more 
or less constant in all three model sequences, although the gas 
pressure changes by almost one order of magnitude. For lower 
[M/H] down to [M/H]s;-3.5, the optical depths start to decrease 
slowly, while the size of the visible dust particles is increasing 
by about one half order of magnitude. This amplitude becomes 
stronger for closer proximities to the convection zone, which in 
our case concerns the low metallicity models of the Teff=2000K 
and log(g)=3.0 sequence. In general, this means larger particles 
are present in the visible dense cloud layers for decreasing metal- 
licities down to [M/H]a!-3.5. The mean particle sizes increase up 
to several 10//m, making it once more difficult to observe the ma- 
jor fraction of the dust in the optically thin parts of the cloud. 

The outward shift of the dust-to-gas ratio maximum below 
[M/H]^-3.5 amplifies the decreasing abundances, which results 
in a strong decline of the visible particle size, while the optical 
depth is continuously decreasing with the metallicity. The opti- 
cal depth rise for increasing metal abundances becomes stronger 
for lower effective temperatures as the amount of dust becomes 
much larger and the dust opacity is increasingly affected by solid 
iron. 

In order to show trends for the amount and average size of 
visible particles, it is best to choose a mean optical depth which 
focuses on those wavelengths at which the dust opacity has a 
strong impact, i.e., where t = 1 is located within the cloud. 
From Fig.|7] one can see that tb is the best applicable optical 
depth of our five examples. Hence, Fig.[8]depicts the dust-to-gas 



ratio and mean grain size at tb = 1 for varied metallicities. From 
Fig.|2]it is obvious that the dust clouds for Teif=2500K are not 
strongly enough developed, which is why r = 1 is only reached 
below the clouds for all considered metallicities. Hence, there is 
neither a dust-to-gas ratio nor a grain size for Teff=2500K sam- 
pled in Fig. [8] For the same reason, there are no data points for 
Teff=2000K, log(g)=3.0 and [M/H]<-4.0 and the Teff=2000K, 
log(g)=5.0 and [M/H]=-6.0 model. Both curves shown for the 
dust-to-gas ratio back the simple assumption that the ratio is de- 
creasing with metallicity. However, the visible dust-to-gas ratio 
in both model sequences drops by merely two orders of magni- 
tude between [M/H]=+1.0 and [M/H]=-3.5 (compare Sec.WM. 
This apparent discrepancy can explained by the cooler atmo- 
spheres for lower metallicities, which enable the dust cloud to 
form in atmosphere layers of higher gas pressure, which involves 
a much more efficient dust formation. This increase in efficiency 
is able to slow down the decrease of the dust-to-gas ratio with 
the metallicity. Only for metallicities below [M/H]=-3.5 the at- 
mospheres will not become significantly cooler for lower metal- 
licities. Thus, the cloud is no longer able to sink any deeper and 
the dust-to-gas ratio is finally dropping fast with metallicity, un- 
til the gas is no longer supersaturated and no dust is able to form. 



In contrast to the visible dust-to-gas ratio, the mean 
grain size at rB = l is increasing for lower metal abundances. 
Comparing Fig.|2] and Fig.|7] one finds that this is due to the 
higher gas densities at low metallicity clouds. In addition, the 
close proximity of the convection zone yields additional growth 
of the mean grain size for the log(g)=3.0 sequence, which is vis- 
ible for [M/H]<-2.5. For a higher surface gravity (log(g)=5.0) 
this is of lesser concern, as it is only visible in the spike at 
[M/H]=-3.5. For lower metallicities, the dust cloud is no longer 
shifted to higher pressures. Therefore, the dust growth becomes 
purely affected by the decreasing abundance and the visible par- 
ticles become smaller and smaller. 
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Fig. 9. Spectral sequences. The red crosses indicate the photometric fluxes in J, H and K bands as they would be observed by 2MASS 
with the respective J, H and K filters. See also Tab. 2. Note that the flux scales difi'er between individual plots. 



5. Results: Spectral appearance 

5.1. Spectral energy distribution for varied metallicities 

The model spectra sequences, including additional spectra for 
[M/H]=-0.5, are shown in Fig.|9] For a better comparison of the 
shape of the spectra, they were binned to a lower resolution. For 
convenience. Tab. 2 contains photometric flux data for the three 
model sequences. We add a note of caution as these results might 
suggest an easy comparison with other atmosphere simulations. 
[HelHng et al. (2008a) did show that different approaches of mod- 
elling dust in substellar atmospheres can yield different results 
already for quantities like mean grain size and dust-to-gas ratio 
(see their Fig. 2) that determine the dust opacity in the radia- 
tive transfer calculations. Their comparison of synthetic colours 
suggests an uncertainty between 10%-30% in the 2MASS JHKs 
colours amongst the different atmosphere simulations. 

The high metallicity models contain strong molecular ab- 
sorption features. Most prominent are the water and carbon- 
monoxide bands between 1.0 and 3.5/im. For [M/H]=-1.5 these 
bands have almost completely vanished in the high log(g) mod- 
els. In the low log(g) models, the same happens for slightly 
lower metallicities. The overall spectra start become much 
bluer, as collision induced absorption (CIA) becomes signifi- 
cant. However, the dust opacity is still strong enough to con- 
ceal atmospheric regions of strong CIA and therefore compen- 



sate its growing influence below 1.0/im and to cause a strong 
absorption between 1 .0 and 2.0/im. Hence, the dust cloud delays 
the trend of the spectra appearing bluer Reduction of the metal 
abundances starts to change this, however In the [M/H]=-2.5 
models, the dust cloud has already become way too transpar- 
ent in the infrared to cover the strong CIA opacity below. The 
impact of CIA is by orders larger than the impact of the dust. 
Despite that, the dust opacity is far from neglegible. Especially 
at wavelengths shorter than 1 .5fj.m the spectrum is still stronly 
affected by the dust. While all other opacities vanish fast with 
the metallicities, the alkali lines in the optical and near-infrared 
are getting stronger down to [M/H]=-4.0. This is because the 
view into deeper atmospheric layers becomes possible, where 
the pressure and density ar e higher, causi ng higher collision rates 
and thereby stronger lines (lJohnasll2007l) . For even lower metal- 
licities, the dust opacity vanishes and the alkali lines get weaker. 
The most extreme metallicity examples ([M/H]=-6.0) are almost 
featureless in their spectra. The dust opacity is finally insignifi- 
cant, while the spectral shape is alone determined by CIA. 



Unfortunately, none of our model spectra show any sig- 
nificantly strong dust opacity features in the near- and mid- 
IR, which might provide a good tool to obs ervationally test 
the dust properties other than by polarimetry dSenguptal 120081; 
iLecavelier Pes Etangs et alj|2008l) . 
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Fig. 5. Dust-to-gas mass ratio pd/Pgas, demonstrating the dust 
content for three model sequences of varying metallicity 
[M/H]. Top: Teff=2000K, log(g)=5.0; Center: Teff=2500K, 
log(g)=5.0; Buttom: Teff=2000K, log(g)=3.0 
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Table 2. Exemplary logarithmic photometric fluxes for JHK- 
bands as seen by 2MASS and UKIRT. 
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Fig. 7. Optical depth t (upper panels) and mean grain size 
(a) (lower panels) at the pressure of the dust-to-gas ra- 
tio maximum f (max{pj/pg)) for varied metallicities. Top: 
Teff=2000K, log(g)=5.0; Center: Teff=2500K, log(g)=5.0; 
Bottom: Teft=2000K, log(g)=3.0 Five optical depths are stud- 
ied: the optical depth at Llfiva (ti 2^in), the Planck-weighted 
mean (tb), weighted means for the mean intensity (rj) and 
surface-normal flux (tf) and the Rosseland-mean (t,oss)- The 
metallicity is displayed at the respective data points. 



5.2. Impact of varied relative abundances 

Furthermore, we have studied the influence of diverging solar el- 
ement abundance data sets. Our standard d ata set for the presen t 
model grid corresponds to the results of Grevesse et al.l ( |1992|) . 
hereafter called Grevesse (old). We compare th e se res ult to mod- 
els, using the newer results bv lGrevesse et a n (I2007h . hereafter 
Grevesse ( new). The third configuration is in accordance with the 
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Fig. 11. Deviation between the respective spectra for the three 
considered solar abundance sets. 



Fig. 8. dust-to-gas ratio (upper panel) and mean grain size (lower 
panels) at tb = 1 as a function of the metallicity. Note that the 
dust cloud is not optically thick for Teff=2500K and, hence, is 
located above tb - I, which is why there are no data points 
for this effective temperature. For the same reason, there are 
no data point for Teff=2000K, log(g)=3.0 and [M/H]<-4.0. and 
Teff=2000K, log(g)=5.0 and [M/H]<-6.0. 
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Fig. 10. Temperature and dust cloud structure for three differ- 
ent solar abundance data sets: Grevesse (old) corresponds to 
iGrevesse et al.l (Il992h. Grevesse (new) t o l&evesse et al] (l2007h 
and Caffau to the IGrevesse et al.1 (l2007l) d ata, but using the re- 



vised O-abundance by Caffau et al. ( 20081) 



Grevesse (new) models, e xcept that we use th e revised oxygen 
abundance determined by ICaffau et al.l (l2008l) . in the following 
shortened to Caffau. 

In general, the differences between our results for the three 
configurations remain extremely small, concer ning the atmo- 
sphere and dust cloud structures. Compared to IGrevesse et al.l 



( I992I) most element abundance values have been adjusted in 
IGrevesse et alj ( l2007h (Grevesse (new)). However, although the 
physics considered for the determination the solar photospheric 
abundances may have been refined, their basis is quite the same. 
Hence, the values of both sets do not differ considerably for the 
majority of elements. The Grevesse (new) configuration con- 
tains a mostly smaller amount of the more abundant heavy el- 
ements. The consequence is a smaller gas opacity, which, in 
turn, results in few Kelvins cooler atmosphere, see uppermost 
panel Fig.[TOl The oxygen abundance of Caffau lies in between 
both Grevesse configurations. Thus, the respective model prop- 
erties are in between as well, though somewhat more resem- 
bling Grevesse (old) with its even higher oxygen abundance. 
Concerning the dust, the slightly reduced molecular abundances 
yield a similar result as a decrease of the metalUcity, but on a 
much smaller scale. The whole dust cloud is marginally shifted 
inwards from Grevesse (old) over Caffau to Grevesse (new). The 
resulting higher gas density at the cloud balances the reduced 
abundances with respect to the concentration of reactants, which 
contribute to the dust growth. Therefore, the mean grain sizes 
and number densities remain almost unaffected. The conclusion 
is that our integrated dust cloud / atmosphere model structures 
are very stable against slight abundance variations. 

The corresponding spectra (Fig.lTTI) can differ remarkably, 
especially in the near-infrared. Hence, in contrast to the dust, 
the photometry of the models depends significantly on the abun- 
dance data set, possibly much stronger than our limited test 
shows. This is particularly interesting fo r studies on the low- 
metallicity class of substellar subdwarfs ( Burgasser et al.ll2009l) 
as the interior element abundance pattern can deviate from the 
solar abundances already for moderate subsolar metallicities of 
[M/H]< -2.0 caused by deviations in the star formation his- 
tory, which has b een observed in e.g., the Carina dwarf galaxy 
jKoch et alj|2008[) . A simple scaling of the element abundances 
might not be sufficient. The interpretation of low-metallicity 
substellar objects, therefore, represents a multi-fold challenge: 
These object are suggested to still produce dust in their atmo- 
spheres, hence, deplete the elements in an inhomogeneous way. 
Additionally, the subdwarfs might exhibit a strongly non-solar 
element abundance pattern depending on the history of their star- 
formation region. 

In order to estimate the influence of extreme abundance pat- 
terns like in the early universe on the atmosphere models, we 
have calculated such a model for Teff=2000K and log(g)=5.0, 
featuring an observed ea r ly uni verse element abundance pattern 
reported by iFrebel et al.l (l2008l) . Therein, most heavy elements 
are considerably less abundant than in the solar environment. On 
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Fig. 12. Comparison of a standard spectrum for Teff =2000K, 
log(g)=5.0, [M/H]=-4.0 (black) and a spectrum for the same ef- 
fective te mperature and surf ace gravity but applying the abun- 
dances of iFrebel et all (|2008|) (red). 

the average, the abundances are comparable to [M/H]<-4.0 for a 
plain scaling of the solar abundances, but the relative abundances 
have changed noticeably. Most important are the elements C, N 
and O which are highly abundant compared to higher order ele- 
ments and, therefore, resemble [M/H]a;-2.0. All this strongly in- 
fluences the chemical equilibrium in the atmosphere. Especially, 
the carbon/oxygen ratio of the given abundance data set has 
changed to a degree that the chemistry is almost tipped to a car- 
bon dominated chemistry ([C/O] ~ 1.0). Much of the carbon is 
already bound in methane because of the increased gas density, 
freeing up a large amount of oxygen. For this reason, strong wa- 
ter bands are present in the near infrared (Fig.fTSI). The amount 
of dust within this atmosphere is of the same order of magni- 
tude as we find it for [M/H]a!-4.0...-4.5 by simple scaling of the 
solar abundances. Merely, some details of the radial dust cloud 
structure are washed out but are still present. However, the early 
universe abundances feature significantly less iron. Hence, the 
mean opacity of the dust grains is strongly reduced as much less 
iron condenses. In general, the overall dust cloud characteristics 
as described in Sect. |4]remain. 

We have shown that the dust cloud is very resistant against 
variations of individual abundances, while the spectra can sensi- 
tively reflect such abundance variations. That is especially true 
for the radical cha nge in the relative abundances reported by 
iFrebel et all (|2008[) which reflect the influence of individual su- 
pernovae in the early universe. The richness of C, N and O yields 
large amounts of CO, H2O, CH4 and N2. Of them, especially the 
H2O leaves strong marks in the near infrared at 1.15, 1.4 and 
l.Sjum. The much weaker dust opacity due to the weaker iron 
contribution permits a deeper view into the atmosphere which, 
combined with slightly higher alkali abundances, results in sig- 
nificantly stronger alkali lines, involving a redistribution of red 
optical flux into the near infrared. In addition, the deeper look 
into the atmosphere provides stronger TiO and FeH features 
shortwards of ljum and slight adjustments to the near-infrared 
CIA opacity. 

6. Conclusions 

The results of this paper on dust formation in low metallicity 
substellar objects are of interest for a variety of objects or scenar- 
ios involving dust at similarely low metalicities, like for instance 



star formation. Our results provide first theoretically obtained in- 
formation on low-metallicity dust characteristics that determine 
the opacity, and we, therefore, generally conclude that 

- Dust can form at very low metallicities but the dust-to-gas 
ratio decreases with decreasing metallicity. 

- The mean grain size decrease with decreasing metalicities 
for a given local temperature and density leaving behind a 
less depleted gas phase then for higher metalicities. 

- The dust material composition changes from Mg/Si-silicate- 
dominated to MgO/Fe-dominated with decreasing metallic- 
ity at about [M/H]=-4.0. 

More specifically, our investigation suggests that dust clouds 
are present throughout the late M- and L-dwarf regime and into 
the regime of hot giant gas planets. Dust is able to form even 
for extremely low metal abundances. Only a combination of 
these low element abundances with higher effective tempera- 
tures or very low surface gravities are able to overturn the strong 
gas phase supersaturation and thereby dissolve the dust cloud. 
Therefore, only the youngest and most massive among brown 
dwarfs and giant gas planets could resist dust formation. 

The dust clouds in our models are shifted inwards for de- 
creasing effective temperature, increasing surface gravity and 
decreasing metallicity. This shift is responsible for the contin- 
uing strong presence of the clouds, as it yields a higher gas 
density, which is in parts able to compensate for the decreas- 
ing metal abundances with [M/H]. Therefore, the maximum dust 
content in our models does not linearly scale with the metallicity, 
but decreases slower than might be expected. 

The convective overshooting becomes more efficient in hot, 
low-gravity objects and also with increasing element abun- 
dances. Hence, young, hot giant planets and brown dwarfs 
are likely to have high-altitude clouds compared to T- or Y- 
dwarfs where the clouds disappear below visibility. Thus, dust 
clouds may turn out as a good tool to distinguish planets 
from substellar companion by a pronounced polarisation sig- 
nal in the planetary case. Indeed the haze-like high-altitude 
clouds of our models support the observational findings of 
Richardson et al. (2007 ) and Pon t etjil. (2 008) w hich were con- 
firmed by iLecavelier Des Etangs et al] (|2008|) and ISenguptal 
(.200&) . based on their polarimetry studies. 

The gas phase at the dust clouds is typically strongly de- 
pleted. In all three model sequences we find a metallicity value 
for which this depletion compared to the deep interior model 
abundances reaches a maximum. For higher and lower metallic- 
ities of the same sequence, the depletion is less strong. Because 
of the strong depletion, which is already present at the upper, 
optically thin cloud, it is risky to determine element abundances 
from absorption lines, as the visible abundances may not reflect 
the deep interior value. 

We have observed a shift of the composition of the dust 
grains. The upper part of the dust cloud, which is dominated by 
silicates, in particular in the case of solar-like metallicities, van- 
ishes due to a developing temperature inversion with decreas- 
ing metallicity. This becomes more severe if the surface gravity 
is getting stronger While dust grains at the bottom of the dust 
cloud are almost fully made of AI2O3 [s] for abundances around 
the solar value, a predominantly Fe[s] volume fraction is found 
in the extremely low metallicity cloud bases, i.e., for the oldest 
amongst the brown dwarfs. High temperatures and surface grav- 
ities amplify this change. 

Dust features are presently not distinguishable in the syn- 
thetic spectra, as the dust opacity features in the near- to mid- 
infrared are too weak and too strongly convolved with the 
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spectral energy distribution of deeper atmosphere layers and 
molecule bands from above the dust cloud. However, the spec- 
tral dust features are extremely reliant on the lattice structure of 
individual dust grains. Unfortunately, no present dust model is 
able to properly treat this issue. Despite the lack of directly ob- 
servable features, the dust clouds remain sufficiently dense to af- 
fect the near- and mid-infrared spectra even below [M/H]=-2.5. 
Adjustments to the code, in order to improve our fitting results, 
like the adaption of a new EOS routine are almost complete. 
Subsequently, a new model grid will be produced and be avail- 
able for pubUc use. 
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